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DEVELOPMENT  OF  A  ROTARY  WING  NAVIER-STOKES  CFD  CODE 

BASED  ON  TLNS3D  CODE 

Hong  Hu 

Hampton  University 
Hampton,  Virginia  23668 


SUMMARY 


This  report  documents  the  development  of  a  multi-grid  Navier-Stokes  code  for  rotary 
wing  calculations,  TLNS3DR.  The  TLNS3DR  is  based  on  the  fixed-wing  TLNS3D  code 
by  addding  a  rotation  term  and  making  necessary  modifications  on  boundary  conditions. 
The  TLNS3D  code  solves  the  thin-layer  Navier-Stokes  equations  using  an  explicit  multi¬ 
stage  Runge-Kutta  type  of  time  stepping  with  multigrid  method.  The  TLNS3DR  code 
re-formulates  the  TLNS3D  code  in  the  rotor  blade-fixed  moving  frame  of  reference  and 
solves  relative  motion  of  fluids. 

In  this  report,  the  formulation  of  the  Navier-Stokes  equations  in  the  moving  frame 
of  reference  in  terms  of  relative  velocity  is  given.  Several  numerical  examples  are  then 
presented.  The  results  of  rotor  blade  in  hover  calculated  by  the  new  TLNS3DR  code 
seems  quantitatively  correct.  The  computation  also  shows  that  this  multigrid  code  is 
efficient. 
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1.  INTRODUCTION 


The  accurate  and  efficient  calculation  of  rotor  flow  fields  poses  a  particularly  chanllege 
problem  in  the  computational  helicopter  aerodynamics  community.  The  accurate  calcula¬ 
tion  of  the  aerodynamic  forces  on  the  rotor  requires  a  correct  prediction  of  the  vortex  wake, 
including  the  blade  trailing  wake  as  well  as  the  interaction  of  the  tip  vortex  of  one  rotor 
blade  with  the  following  rotor  blade.  The  generation  of  the  free-wake  and  its  interaction 
with  the  blade  make  the  helicopter  aerodynamic  computation  extremely  complex. 

The  computational  methods  for  calculating  helicopter  rotor  flow  fields  can  be  divided 
into  two  types:  integral  equation  methods  (panel  methods)  and  finite-difference/finite- 
volume  methods.  Integral  equation  methods  for  both  incompressible  and  compressible 
rotor  flows  are  developed  by  several  investigators1-5.  The  methods  are  able  to  treat 
rotor  blade  and  vortex  wake  accurately  and  efficiently.  However,  they  are  restricted  to 
linear,  low- speed,  subsonic,  or  at  most,  to  nonlinear  flows  with  weak  shocks.  This  is 
a  very  serious  limitation,  because  modern  rotor  blades  operate  in  the  transonic  regime, 
where  the  nonlinear  compressibility  is  important  and  the  strong  shocks  can  exist  in  the 
flow  field.  Finite-difference/finite-volume  methods  based  on  the  full-potential6’7,  Euler8-10 
and  Navier-Stokes9,10  (NS)  equations  with  a  single  grid  for  rotor  flow  fields  are  developed 
during  the  past  few  years,  while  more  recently  the  calculations  are  mainly  concerned  with 
various  NS  equations.  These  methods,  particularly  NS  including  thin-layer  NS  methods, 
can  treat  compressible  flows  with  shocks,  including  strong  shocks,  in  the  region  of  the 
rotor  blade  accurately.  The  effects  of  free-wake  are  captured  fairly  accurately  as  a  part 
of  the  overall  flowfield.  However,  the  calculation  of  NS  equations  on  a  fine,  single  grid  is 
expensive. 

On  the  other  hand,  a  multigrid  NS  Computational  Fluid  Dynamics  (CFD)  computer 
code,  TLNS3D,  is  developed11  for  the  fixed-wing  motion.  The  TLNS3D  code  solves  the 
thin-layer  Navier-Stokes  equations  using  an  explicit  multistage  Runge-Kutta  type  of  time¬ 
stepping  with  multigrid  method.  The  philosophy  of  multigrid  method  is  the  use  of  suc¬ 
cessively  coarser  grids  to  compute  corrections  to  a  fine  grid  solution.  The  use  of  multigrid 
strategy  results  in  a  much  faster  convergence  than  single-grid  method.  The  efficiency  of  the 
multigrid  method  using  five-stage  Runge-Kutta  time-stepping  is  demonstrated11  through 
computations  for  transonic  flows  over  fixed  wings.  The  code  is  shown  to  result  in  a  sig¬ 
nificant  reduction  in  CPU  time  as  compared  to  the  single-grid  multistage  time-stepping 
method. 

The  purpose  of  the  present  work  is  to  develop  a  rotary-wing  multigrid  Navier-Stokes 
code,  TLNS3DR,  by  modifying  the  TLNS3D  code.  The  new  code  re-formulates  the 
TLNS3D  code  in  the  rotor  blade-fixed  moving  frame  of  reference  and  solves  relative  mo¬ 
tion  of  fluids.  In  this  report,  the  formulation  of  the  Navier-Stokes  equations  in  the  moving 
frame  of  reference  in  terms  of  relative  velocity  is  given.  Several  numerical  examples  are 
then  presented.  The  results  of  rotor  blade  in  hover  calculated  by  the  new  TLNS3DR 
code  seems  quantitatively  correct.  The  computation  also  shows  that  this  multigrid  code 
is  efficient. 


2.  THE  NAVIER-STOKES  EQUATIONS 
2.1  NS  Equations  in  Space-fixed  Frame  of  Reference 
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The  basic  equations  under  consideration  are  unsteady  Navier-Stokes  equations.  In  the 
space-fixed  frame  of  reference,  the  dimensional  form  of  the  NS  equations  in  the  conservative 
form  is  given  by 

%  +  V*  ■  (p'V')  =  0  (1) 

Jr(/>*r*)  +  v*  ■  +  P*I)  -  v*  •  r  =  o  (2) 

r\ 

— (/e*)  +  V*  ■  (p*h*V*)  -  V*  •  (f*  ■  V*  -  f)  =  0  (3) 


with 


P*  =  p'{ 7  -  l)(c* 


y*2 


2 


) 


(4) 


h* 


IP* 

(7  -  1  )p* 


(5) 


where  p*  is  the  density,  V*  the  velocity  vector,  t*  the  time,  p*  the  pressure,  I  the  identity 
matrix  tensor,  f  the  viscous  shear  stress  tensor,  e*  the  total  energy  given  by  e*  =  e*+V*2/2 
with  e*  as  internal  energy,  h*  the  total  enthalpy,  and  q*  the  heat  conduction  given  by 

q*  =  — &V*T*  where  T*  is  the  temperature.  The  superscript,  *,  denotes  dimensional 
quantities. 


Using  the  free-stream  density  (p^),  the  free-stream  speed  of  sound  divided  by  the 
square-root  of  the  gas  specific  heat  ratio  (a^/^/y),  and  wing  root-  chord  length  (c),  the 
above  equations  take  the  following  form: 


^+V.<J>V)  =  0 


(6) 


with 


ps 

~{PV)  +  V  •  (pVV  +  pi)  -  V  •  f  =  0 

A 

-(pe)  +  V-(phV)-V-(f-V-?)  =  0 


V2 

p  =  p(7-l)(e  -  — ) 


h  = 


7  p  V2 

(7-1  )p  +  ~2~ 


(7) 

(8) 

(9) 

(10) 


where  p  =  p* / V  =  cV*;  V  =  V*s/y/ 

T  =  f*/p*oo ;  e  =  e*7/«*L;  h  =  h^/a*2^, 

refers  to  free-stream  condition. 


<0;  t  =  faSoAVr);  P  =  P*7/p*oc^*lo  =  P*/Pt o; 
9  =  and  the  subscript,  00, 
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Equations  (6),  (7)  and  (8)  along  with  Eqs.  (9)  and  (10)  are  Navier-Stokes  equations 
in  non-dimensional  form,  where  f  and  q  are  given  by 


and 


R, 


9  —  (Qxi  Qy: 9z) 


&uj  \  2  du  £  -j 
dx(  3  t':' dxk 

(ii) 

RePr(7  ~  1) 

(12) 

respectively,  where  fx  is  the  non-dimensional  viscous  coefficient  given  by  fi  =  fx*//x ^  = 
T3/2(  1  +  C)/(T  +  C),  and  C  =  0.4317;  T  is  the  non-dimensional  temperature  given  by 
T  =  T*/T*\  Re  is  the  Reynolds  number  given  by  Re  =  is  the  free- 

stream  Macb  number  given  by  =  V^Ja*^  and  PT  is  the  Prandtl  number  given  by 
Pr  =  Cpfi^/k^  with  cp  as  gas  specfic  heat  and  k*^  as  free-stream  conductivity. 


2.2  NS  Equations  in  Blade-fixed  Frame  of  Reference 


For  a  general  motion  of  a  rotor  blade,  the  governing  equations  are  simple  to  solve 
if  a  rotor  blade-fixed  moving  frame  of  reference  formulation  is  used.  In  addition  to  the 
space-fixed  frame  of  reference  OXYZ,  a  blade-fixed  moving  frame  of  reference  oxyz  is 
introduced  as  shown  in  Figure  1.  The  moving  frame  of  reference  oxyz  is  translating  at 

a  velocity  of  V0(t )  and  rotating  around  a  pivot  point  (rotor  blade  axis:  rp(xp,  yp,  zp  j) 
at  an  angular  velocity  of  f2(i).  The  relation  for  the  non-dimensional  absolute  velocity 
( V  =  (w,u,iu)),  relative  velocity  ( Vr  =  (ur,vr,wr))  and  transformation  velocity  ( Vt  = 
(ut,  vt ,  wt )  =  V0  +  Ye  —  Vo  +  x  r)  is  given  by 


V  =  Vr  +  V0  +  fixr  (13) 

where  r  is  the  position  vector  measured  in  the  blade-fixed  moving  frame  of  reference  relative 

to  rp.  It  should  be  noted  that  the  transformation  velocity  (Vt)  includes  both  translation 
— * 

(V0)  and  rotation  (fi  x  f)  velocities  and  therefore  this  formulation  is  general  for  both 
bovering  and  forwarding  motion  of  the  helicopter. 


To  express  Eqs.  (6),  (7)  and  (8)  in  terms  of  the  blade-fixed  moving  frame  of  reference, 
the  relations  of  substantial  and  local  derivatives  are  used: 


Db  D'b 

Dt  ~  Dt' 

(14) 

9b  d'b  ^ 

dt~  dt'  Vt ' V 

(15) 

DB  D'B  - 

Dt  =  Dt' 

(16) 

dB 

dt 

d'B  - 

=  ^--Vt-XB  +  QxB 
dt' 

(U) 
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where  the  prime  (')  denotes  the  operation  with  respect  to  the  blade-  fixed  moving  frame 

— * 

of  reference;  b  is  a  scalar  and  B  is  a  vector. 

By  using  Eqs.  (13)  through  (17),  Eqs.(6),  (7)  and  (8)  become 


+  V  •  (pVr)  =  0  (18) 

^i(pVr)  +  V  •  (pVrVr  +  pi)  -  V  •  fr  =  -pSt  (19) 

^(per)  +  V  •  (phrVr)  ~  V  •  (fr  •  Vr  -  q) 

=  -p[Vr  •  a0  +  (fi  x  r)  ■  a0  +  V0  ■  (at  -  0  x  Vr)  (20) 

+  x  r)  +  x  ^  •  (^r  x  0] 

with 


V2  V? 

p=(7-lWer-^-  +  -f) 

(21) 

,  7 p  .  n  t? 

r  pin  —  i)  '  2  2 

(22) 

V  N  DV  D'Vr 

at  —  (°tx>  atyi  atz)  — 

d'O  -  _ 

=  a0  +  {—xr)  +  2(nxVr)  +  nx(Qxr) 

(23) 

=  Mooy/jp  wduri  i  duTj  N  2  c  durk] 

Re  dxj  dxi  3  t,J  dxk 

(24) 

where  a0  is  the  acceleration  of  the  translation  velocity  given  by  a0  =  DV0/Dt  and  hr  is  the 
total  rothalpy.  Equations  (18),  (19)  and  (20)  along  with  Eqs.  (21)  through  (24)  are  set  of 
Navier-Stokes  equations  in  the  moving  frame  of  reference  in  terms  of  relative  velocity. 


In  the  cartesian  coordinates,  Eqs.  (18),  (19)  and  (20)  are  rewritten  in  the  matrix  form 


as: 


dl'  +  fjx<F'  ^"r>  +  dy1'^’  6"r>  +  dzHr  ^ 


(25) 


Ur  =  [p,pur,pvr,pwr,perY 
Fr  =  [pur,  pu2r  -f  p,  purvr,  purwr,  purhrY 
Gr  =  [pvr,  pvrur,  pv2T  +  p,  pvrwr,  pVThrY 


with 
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(26) 

(27) 

(28) 


Hr  =  [pwr,  pWrUT ,  pWrVr,  pw2r  +  p,  pwrhrY 

(29) 

Fvr  —  [0,  Trxx  ,  Trxy  ,  Trxz  ,  UrTrxx  +  VrTrxy  4  wrTrxz  ~  Px) 

(30) 

Gvr  =  [0,  Tryxi  Tryy  i  Tryz  i  urTryx  4"  vrTryy  4"  wrTryz  ~  Qy] 

(31) 

Hvr  —  [0,  Trzx  ,  T~rzy  i  Trzz  i  urTrzx  4  vrTrzy  4"  WrTTzz  ~  <Jz\ 

(32) 

S  =  [0,  -patx ,  -paty ,  -patz , 

-  p[Vr  ■  a0  +  (O  x  r)  ■  a0  +  V0  •  (at  -  Cl  x  Vr) 

(33) 

+  Vr'(~dF  x  ^)  +  (Q  x  r)-(-^r  x  *01] 


where  Ur  is  the  field  vector;  Fr,  Gr  and  Hr  are  inviscid  flux  vectors;  Fvr,  Gvr  and  Hvr 
are  viscous  flux  vectors;  and  S  is  the  source  vector  due  to  the  motion  of  the  moving  frame 
of  reference.  All  field  and  flux  vectors  and  viscous  shear  stresses  are  based  on  the  relative 
velocity  field.  The  subscript,  r,  denotes  relative  quantity. 

2.3  NS  Equations  in  Computational  Domain 


The  Navier-Stokes  equations  given  by  Eq.  (25)  are  then  specialized  to  a  body  fit¬ 
ted  coordinate  system  (£,?/,()  in  computational  domain,  where  £,  p  and  (  represent  the 
streamwise,  normal  ana  spanwise  directions,  respectively.  The  thin-layer  assumption  is 
employed  by  keeping  viscous  diffusion  term  in  the  normal  (p)  direction  only,  since  the 
dominant  viscous  effects  arise  from  viscous  diffusion  normal  to  the  body  surface  for  high- 
Reynolds-number  turbulent  flows.  For  the  time-independent  (£,  p,  ()  coordinate  system  in 
the  blade-fixed  moving  frame  of  reference,  Eq.  (25)  becomes 


§(J-‘^  +  aF' 


r  dGr  dHr  _  dGvr 

d£  dp  ^  dC,  dr] 


+  J-'S 


with 

Fr  =  J  1(£xFr  +  £yGr  +  £zHr) 

=  J~1[pur,  purur  +  £xp,  purvr  +  Zyp,purwr  +  £zp,purhrY 
Gr  =  J  1(pxFr  +  rjyGr  +  pzHr) 

=  J-1  [pvr,  pvrur  +  T]xp,  pvrvr  +  T]yp,  pVrWr  +  T]zp,  pVrhrY 

Hr  =  J  x(Ci +  (yGr  +  (zHr) 

=  J-1  [pwr,  pwrur  +  (xp,  pwrvr  +  CyPi  pwrwr  +  C zp ,  pwrhrY 


Gyr  *7  (p XFVr  PyG vr  ”1”  PzHyfa 

-  J  - ^ fa  Urv  +  Pxfai  favrv  +  Py<f>2 ,  fawrr,  +  Pz<j>  2,  +  fafa  ] 


where 


Ur  —  (iWr  “I-  ^yVr  “I" 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 
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Vr  =  rjxur  +  r/yvr  +  T]zwr  (40) 

Wr  =  Cx^r  4"  Cy^r  4~  Cz^r  (41) 

fa  =  vl  +  r}2y  +  V2z  (42) 

fa  =  ^(VxUrv  +  T]yVrjl  +  TjzWr7])  (43) 

q2  =  u2r  +  v  2  -f  w2r  (45) 


and  J-1  is  the  inverse  Jacobian  of  the  transformation.  It  should  be  noted  that  the  viscosity 
(fi)  and  conductivity  ( k )  are  replaced  by  fie  and  ke,  respectively,  to  include  the  effect  of 
turbulence,  where  e  is  turbulence  edge  viscosity  and  e  is  tuebulence  edge  conductivity.  The 
turbulence  model  used  here  is  the  algebraic  model  of  Baldwin  and  Lomax12. 

Equation  (34)  is  the  most  general  form  of  thin-layer  NS  equations  in  the  blade-fixed 
moving  frame  of  reference  in  terms  of  relative  velocity.  The  5-term  in  Eq.  (34)  is  the 

source  term  contributed  from  the  time-dependent  rotation  (0(f))  and  translation  ( V0{t )) 
of  the  rotor  blade.  For  hovering  motion  which  is  the  case  considered  here,  and  where 

V0  =  0  and  0  =  0 i  —  0 j  +  Ok  with  0  =  constant,  the  5-term  in  Eq.  (34),  or  Eq.  (33), 
reduces  to 

5  =  p[ 0,  (x  —  xp) O2  -f  20ior,  0,  (z  —  zp) 02  —  20ur,  0]1  (46) 

where  (xp,yp,zp)  is  the  blade  rotation  axis  as  mentioned  earlier. 


3.  NUMERICAL  SCHEME 


The  present  work  is  based  on  the  fixed-wing  TLNS3D  code11.  The  difference  of  the 
present  NS  formulation  from  the  fixed-wing  NS  formulation  of  TLNS3D  code  is  due  to 
the  source  term,  5.  Hence  the  numerical  scheme  follows  exactly  the  one  of  the  TLNS3D 

code  with  the  addition  of  the  source  term  (5)  into  the  code,  along  with  the  necessary 
modification  of  the  boundary  conditions. 

3.1  Finite- Volume  Discretization 


Replacing  the  spatial  derivatives  by  central  differences  in  Eq.  (34),  the  following 
semi-discrete  form  is  obtained: 


Ur)i,j,k  +  {Fri+\/2,j,k  Fri-l/2,j,k ) 

4"  {Gri,j+l/2,k  ~  Gri,j—\/2,k) 

4"  {Hfi,j,k+ 1/2  —  Hri  j  k- 1/2) 

=  (Gvri,j+l/2,k  ~  Gvri,j-l/2,k)  + 


(47) 
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where 


Fri±\/2,j,k  2  ^ri>J>k  ^ri±l,j,k) 

Gri,j±  1/2, k  —  2^ri’j<k  Gri,j±l,k) 

Hri,j,k±l/2  =  +  Hri,j,k±l) 


(48) 

(49) 

(50) 


and  A£,  A rj  and  AC  are  taken  to  be  one.  The  cell  centroidal  point  is  ( i,j,k ),  and  the 
halh-integer  subscripts  refer  to  cell  interface.  The  viscous  term  (Gvrij± i,*)  is  evaluated 
according  to  TLNS3D  code  using  relative  velocity.  Artificial  dissipation  of  TLNS3D  code 
is  used  where  the  velocity  is  repleced  by  the  relative  velocity. 


Using  the  notation  Qr  to  represent  all  inviscid  flux  terms,  Qvr  to  represent  the  viscous 
flux  terms  and  Dr  to  represent  the  artificial  dissipative  fluxes  for  convenience,  Eq.  (47)  is 
rewritten  as 


^(J  Ur)i,j,k  +  Qr  —  Qvr  +  Dr  + 


(51) 


where 


Qr  —  (Fri+1,2,j,k  ~  Fri-l/2,j,k) 


and 


+  (Gri,j+i/2,k  -  Ori,j-l/2,k)  (52) 

+  (Hri,j,k+ 1/2  -  Hrij  k- 1/2) 

Qvr  =  (GVrij+i/ 2,k  ~  ^uri,j-l/2,jt)  (53) 


3.2  Multi-stage  Time-Stepping  Scheme 

The  source  term  ( S )  is  added  into  the  five-stage  Runge-Kutta  time-  stepping  scheme 
used  in  the  TLNS3D  code.  This  time-stepping  scheme  is  described  as  follows: 

U(r0)  = 

U?1  =  I7<0>  -  a,  A;-! Q<°)  +  M'S  -  7iD<°>  -  J~1S] 

U<2>  =  t/<°>  -  02  +  hdfS  ~  7 140)  -  72B'1'  -  J-‘5) 

P<3>  =  U<°>  -  o3  J^-[Q<2)  +  M'S  -  7i D<°>  -  72 D'S  -  73fl<2)  -  /->5) 

u'S  =  u<°>  -  04  fS[Qm  +  PiQm  _  llDm  _  _  l3Dm  _  yiDm  _  7-i  jj 

p(5>  =  (/(o)  _  o5 +fta<-r)  -  7ib<»)  -  12dS  -  73 dS  -  74 i><3)  -  75 dS 

-  j-'s] 

lf(n+l)  =  fj{  5) 


(54) 
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where  the  superscripts,  (n)  and  (n  +  1),  refer  to  time  levels;  and  the  superscripts,  (0) 
through  (5),  refer  to  multistage  level  within  each  time  step.  In  addition  to  the  five- 
stage  time-stepping  scheme,  the  local  time-stepping  and  residual  smoothing  technigues 
are  applied  to  accelerate  the  convergence  of  the  scheme  to  steady  state  for  rotor  hovering 
motion.  The  details  of  the  residual  smoothing  is  discussed  in  Ref.  [11]. 

3.3  Multigrid  Acceleration  Technique 

The  multigrid  scheme  of  TLNS3D  is  employed  for  hovering  motion  calculations.  The 
idea  of  multigrid  is  the  use  of  successively  coarser  grids  to  compute  corrections  to  a  fine 
grid  solution.  The  use  of  multigrid  technique  can  significantly  accelarate  the  convergence 
rate.  In  order  to  provide  optimum  damping  of  the  high  frequency  errors,  the  coefficients 
of  ai  =  1/4,02  =  1/6,03  =  3/8,04  =  1/2,05  =  1,  0j  =  1,02  =  0,03  =  1,04  =  0,05  =  1, 
and  71  =  1,72  =  0,73  =  0.56, 74  =  0,75  =  0.44  are  used  according  to  TLNS3D.  The 
coefficients  on  the  residual  smoothing  are  also  carefully  selected.  The  results  presented  in 
this  report  are  obtained  using  a  V-cycle  multigrid  scheme. 

3.4  Boundary  Conditions  and  Grid 

A  C-0  type  grid  is  employed.  The  grid  is  generated  using  the  WTCO11  code.  Figure 
2  is  the  partial  view  of  C-0  grid  for  the  half-span  rectangle  rotor  blade.  A  mesh  consisting 
of  129  x  65  x  65  points  in  the  streamwise,  spanwise  and  normal  directions,  respectively,  is 
used. 

On  the  blade  surface,  no-slip  and  no-penetration  conditions  are  used  by  setting  relative 

velocity  (Vr)  to  zero.  The  adiabatic  wall  condition  and  zero-normal  pressure  gradient 
condition  at  the  wall  are  also  applied  at  the  blade  surface.  The  farfield  boundaries  are 
treated  by  using  Riemann  invariants  condition  used  in  TLNS3D  code11.  At  the  inboard- 
plane,  near  the  rotation  axis  of  the  blade,  the  condition  of  urjwT  —  ut/wt  =  —(z  —  zp)/(x- 
xp )  is  used  to  guarantee  that  the  radial  velocity  is  zero  at  this  boundary. 

4.  NUMERICAL  EXAMPLES 

Calculations  are  performed  for  a  straight,  untapered  and  untwisted  rectangle  blade 
with  constant  NACA0012  section.  A  single  isolated  blade  is  considered.  The  blade  length 
is  12,  while  the  chord  length  is  1.  Flows  at  three  different  conditions  for  hovering  motion 
are  calculated:  non-lifting  shock-free  subsonic  flow,  non-lifting  transonic  flow  and  lifting 
transonic  flow.  For  the  non-lifting  flows,  the  calculations  for  fixed-wing  motion  under 
equivalent  condition  are  also  made  for  serving  as  reference  solutions.  Additional  extensive 
calculations  are  made  on  the  BERP  blade,  the  UH-60A  blade  and  a  equivalent  rectangle 
blade,  and  results  are  presented  in  Ref.  [13]. 

The  first  example  is  for  non-lifting  shock-free  subsonic  flow.  Figures  3a  and  3b  are  the 
surface  pressure  and  local  Mach  contours,  respectively,  for  fixed-wing  flow  at  M =  0.628 
and  Re  =  3.7  x  106.  All  Mach  contours  presented  here  are  calculated  at  j  =  33  which 
is  about  just  outside  the  boundary  layer.  Figures  3c  and  3d  are  the  surface  and  local 
Mach  contoures,  respectively,  for  hovering  motion  at  Mtip  =  0.628  and  Re  =  3.7  x  106. 
By  comparing  the  results  of  hovering  motion  with  those  of  fixed-wing  motion,  it  is  seen 
that  effects  of  rotation  term,  which  is  J-15  in  NS  equations,  are  quantitatively  correctly 
calculated.  The  solution  for  hovering  motion  looks  reasonably  correct.  Figures  3e  and  3f 
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are  convergence  history  in  terms  of  residual  for  fixed-wing  motion  and  hovering  motion, 
respectively.  A  total  of  400  equivalent  fine  grid  time  steps  (beased  on  the  equivalent  CPU 
time  required)  are  used  in  both  cases.  The  CPU  time  on  the  NASA-LaRC’s  CRAY-YMP 
with  a  single  processor  for  all  cases  presented  here  is  around  6100  seconds,  which  is  about 
0.011  seconds  per  fine  grid  point,  or  27.98  x  10-6  seconds  per  fine  grid  point  per  time  step. 

The  second  example  is  for  non-lifting  transonic  flow.  Figures  4a  and  4b  are  the  surface 
pressure  and  local  Mach  contours,  respectively,  for  fixed-wing  flow  at  M0 0  =  0.84  and 
Re  =  4.97  x  106.  Figures  4c  and  4d  are  the  surface  and  local  Mach  contoures,  respectively, 
for  hovering  motion  at  Mtip  =  0.84  and  Re  —  4.97  x  106.  By  comparing  the  results  of 
hovering  motion  with  those  of  fixed-wing  motion,  it  is  seen  that  effects  of  rotation  term  are 
quantitatively  correctly  calculated  for  this  transonic  flow  case.  The  solution  for  hovering 
motion  looks  reasonably  correct,  where  shock  occurs  in  the  tip  region  of  the  blade.  Figures 
4e  and  4f  are  convergence  history  in  terms  of  residual  for  fixed-wing  motion  and  hovering 
motion,  respectively.  These  two  figures  tell  us  that  the  solution  for  hovering  motion  does 
not  converge  as  fast  as  one  for  fixed-wing  motion,  and  this  fact  may  due  to  the  unsteadiness 
of  the  hovering  motion. 

The  third  numerical  example  is  for  lifting  transonic  flow.  Figures  5a  and  5b  are  the 
surface  and  local  Mach  contoures,  respectively,  for  hovering  motion  at  angle-of-attack  of 
4.2°,  Mtip  =  0.877  and  Re  =  3.93  x  106.  The  solution  looks  reasonably  correct,  where 
shock  occurs  in  the  tip  region  of  the  blade  with  a  reasonable  shock  strength  and  location. 
Figure  5c  is  the  tip  cross  flow  vector  field  where  the  relative  velocity  is  plotted.  From  this 
figure,  it  is  seen  that  the  tip  roll-up  vortex  is  captured  with  the  current  (129  x  65  x  65) 
mesh  size,  and  the  solution  of  the  tip  vortex  looks  quantitatively  correct.  Figures  5d  is 
convergence  history  in  terms  of  residual. 

5.  CONCLUDING  REMARKS 

The  TLNS3D  code  is  modified  into  TLNS3DR  code  for  rotary-wing  calculations.  The 
solutions  for  helicopter  rotor  blade  in  hover  are  obtained  using  the  new  TLNS3DR  code 
and  the  solutions,  including  shock  and  tip  roll-up  vortex,  looks  quantitatively  correct.  The 
accurancy  of  the  new  code  is  going  to  be  tested,  and  the  fine-tune  of  the  code  may  be 
necessary  to  improve  quality  of  the  code.  The  application  of  multigroid  method  in  the 
present  rotary-wing  calculations  shows  that  the  method  is  very  efficient  compared  with 
the  single  grid  method,  although  the  solution  for  hovering  motion  converges  slower  than 
that  for  fixed-wing  motion.  The  code’s  capability  of  capturing  free  tip  vortex  is  promising 
as  a  analysis  tool  in  rotor  tip  design. 
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NACA001 2  Section,  AR=6, 6=0° 

Partial  View  of  C-0  Mesh  for  Rectangle  Blade,  129x65x65  Mesh 


Figure  2.  C-0  mesh. 
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Upper  Surface 


Lower  Surface 


NACA001 2  Section,  Non-Rotating  Rectangle  Blade 
6=a=0°,  Mx=.628,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 


Level  NCP 
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F  -0.299 

E  -0.351 

D  -0.404 

C  -0.457 
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8  -0.667 

7  -0.720 

6  -0.772 

5  -0.825 

4  -0.878 

3  -0.930 

2  -0.983 

1  -1.035 


Figure  3a.  Non-rotating  subsonic  flow,  pressure  coefficient. 
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Upper  Surface 


Lower  Surface 


NACAOOI 2  Section,  Non-Rotating  Rectangle  Blade 


atj=33 

Level  Mach 

U  0.773 

T  0.757 

5  0.741 

R  0.725 
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8  0.421 

7  0.405 

6  0.389 

5  0.373 

4  0.357 

3  0.341 

2  0.325 

1  0.309 


0=a=O°,  M^.628,  Re=3.7x106, 129x65x65  Mesh 


Calculated  Local  Relative  Mach  Contours 


Figure  3b.  Non-rotating  subsonic  flow,  local  Mach  number. 
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Level  NCP 


U  0.704 

T  0.646 

5  0.589 

R  0.531 

Q  0.474 

P  0.416 

O  0.359 

N  0.302 

M  0.244 

L  0.187 

K  0.129 

J  0.072 

1  0.015 
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G  -0.100 

F  -0.158 

E  -0.215 

D  -0.273 

C  -0.330 

B  -0.387 

A  -0.445 

9  -0.502 

8  -0.560 

7  -0.617 

6  -0.674 

5  -0.732 

4  -0.789 

3  -0.847 

2  -0.904 

1  -0.962 


NACA001 2  Section,  Rectangle  Blade  in  Hover 
0=a=O°,  Mtjp=.628,  Q=0.0619,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 


Figure  3c.  Subsonic  flow  in  hover,  pressure  coefficient. 
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at  j =33 

Level  Mach 

U  0.830 
T  0.811 

5  0.791 
R  0.772 
Q  0.752 
P  0.732 
O  0.713 

Upper  Surface  n  0.693 

rr  M  0.674 

L  0.654 
K  0.634 
J  0.615 

1  0.595 

H  0.576 
G  0.556 
F  0.537 
E  0.517 
D  0.497 
C  0.478 
B  0.458 
A  0.439 
9  0.419 

8  0.400 

7  0.380 

6  0.360 

5  0.341 

4  0.321 

3  0.302 

2  0.282 

NACA001 2  Section,  Rectangle  Blade  in  Hover  1  0  263 
0=a=O°,  Mtip=.628,  Q=0.0619,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Local  Relative  Mach  Contours 


Figure  3d.  Subsonic  flow  in  hover,  local  Mach  number. 
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NACA001 2  Section,  Non-Rotating  Rectangle  Blade 
0=a=O°,  M  =.628,  R  =3.7x1 06, 129x65x65  Mesh 

9  oo  8 

Convergence  History  in  terms  of  Residual 


Figure  3e.  Non-rotating  subsonic  flow,  residual. 
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200 


400 

Work 

NACA0012  Section,  Rectangle  Blade  in  Hover 

e=a=0°,  Mtip=.628,  Q=0.0619,  Re=3.7x106, 129x65x65  Mesh 
Convergence  History  in  terms  of  Residual 


Figure  3f.  Subsonic  flow  in  hover,  residual. 
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Level  NCP 


Upper  Surface 


Lower  Surface 


U  0.840 

T  0.772 

5  0.705 

R  0.637 

Q  0.570 

P  0.502 

O  0.435 

N  0.367 
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L  0.232 

K  0.165 

J  0.097 

1  0.030 

H  -0.038 
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D  -0.308 

C  -0.375 
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8  -0.645 

7  -0.713 

6  -0.780 

5  -0.848 

4  -0.915 

3  -0.983 

2  -1.050 
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NACA001 2  Section,  Non-Rotating  Rectangle  Blade 
0=a=O°,  M^.84,  Re=4.97x1 06, 129x65x65  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 


Figure  4a.  Non-rotating  transonic  flow,  pressure  coefficient. 
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Upper  Surface 


Lower  Surface 


NACAOOI 2  Section,  Non-Rotating  Rectangle  Blade 


atj=33 


Level  Mach 
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D  0.753 

C  0.721 

B  0.689 

A  0.657 

9  0.625 

8  0.593 

7  0.562 

6  0.530 

5  0.498 

4  0.466 

3  0.434 

2  0.402 

1  0.370 
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Figure  4b.  Non-rotating  transonic  flow,  local  Mach  number. 
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NACA001 2  Section,  Rectangle  Blade  in  Hover 
0=a=O°,  Mtip=.84,  Q=0.0828,  Re=4.97x106, 129x65x65  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 


Figure  4c.  Transonic  flow  in  hover,  pressure  coefficient. 
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NACA001 2  Section,  Rectangle  Blade  in  Hover 
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Figure  4d.  Transonic  flow  in  hover,  local  Mach  number. 
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NACA001 2  Section,  Non-Rotating  Rectangle  Blade 
0=a=O°,  Mk=. 84,  Re=4.97x106, 129x65x65  Mesh 
Convergence  History  in  terms  of  Residual 


Figure  4e.  Non-rotating  transonic  flow,  residual. 
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Convergence  History  in  terms  of  Residual 


Figure  4f.  Transonic  flow  in  hover,  residual. 
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NACA001 2  Section,  Rectangle  Blade  in  Hover 
0=ct=4.2°)  Mtip=.877,  Q=0.0865,  Re=3.93x106, 129x65x65  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 


Figure  5a.  Lifting  transonic  flow  in  hover,  pressure  coefficient. 
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Figure  5b.  Lifting  transonic  flow  in  hover,  local  Mach  number. 
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NACA001 2  Section,  Rectangle  Blade  in  Hover 
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Figure  5c.  Lifting  transonic  flow  in  hover,  tip  flow  field  investigation. 
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NACA0012  Section,  Rectangle  Blade  in  Hover 
a=4.2°,  Mtip=.877,  Q=0.0865,  Re=3.93x106, 129x65x65  Mesh 
Convergence  History  in  terms  of  Residual 


Figure  5d.  Lifting  transonic  flow  in  hover,  residual. 


30 


LIST  OF  PUBLICATIONS 


None. 


31 


LIST  OF  PARTICIPATING  SCIENTIFIC  PERSONNEL 


1.  Hong  Hu,  Principal  Investigator,  Associate  Professor,  Department  of  Mathematics,  Hamp¬ 
ton  University. 

2.  Terry  G.  Logan,  Graduate  Student,  Department  of  Mathematics,  Hampton  University, 
Degree  Earned:  Master  of  Science  in  Applied  Mathematics,  (the  M.S.  Thesis  was  supported 
by  the  NASA  Grant  NAG-1-1170). 

3.  George  H.C.  Page,  Jr.,  Graduate  Student,  Department  of  Mathematics,  Hampton  Univer¬ 
sity. 


32 


LIST  OF  INVENTIONS 


None. 


33 


